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^ ! Abstract 

We study the effect of finite temperatures and Fermion density 
Q\ ■ on the effective pseudo-scalar-photon vertex induced by the triangle 

diagram. The manifestly covariant calculations show that when the 
Qh| pseudo-scalar mass is much less than the temperature, then there is a 

q; large enhancement of the decay rate. Alternatively, when the temper- 

(D . ature is much higher than the Fermion mass, or the Fermion chemical 

potential is large, the lifetime is enhanced. Other related processes, 
and applications of these results to cosmology and astrophysics, are 
S^ . discussed. 
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1 Introduction 

Goldstone Bosons arise in field theories when some global symmetry is spon- 
taneously broken. A solution of the strong CP problem involves the breaking 
of Peccei-Quinn symmetry and gives rise to a pseudo-scalar called the axion 
(A) [|H], which becomes massive through mixing with the pion. Axions have 
couplings to charged Fermions and are candidate dark matter which could 
play an important role in cosmology . The physics of axions involves two 
intrinsic mass scales — the axion mass, M, and the Fermion mass, m. If 
axions participate in physics inside a heat bath with temperature, T, larger 
than any one of (or both) these scales, then we may expect interesting and 
non-trivial thermal effects. Since the mass of the axion can be much less than 
the cosmic microwave background temperature (T& = 2.7 K), even the cold- 
est environment in the universe may be thermal in this sense. The purpose 
of this paper is to identify the main sources of non-trivial thermal effects. 

We compute the decay rate of the axion into two photons through the 
triangle diagram and find two sources of thermal effects. Stimulated emission 
of final state photons enhances the decay width of the pseudo-scalar, and 
plays the dominant role in the physics of axions for M < T&. Pauli blocking 
of on-shell parts of the Fermions in the loop reduces the decay width when 
T > m, and may have interesting consequences in the physics of supernova?. 
These changes in the effective A77 vertices also contribute to the Primakoff 
effect, 7e — > Ae through a t-channel photon exchange. Thermal effects can 
also boost the multi-photon channels to a very large extent and may require 
a resummation of this whole class of processes. 

One bit of physics is worth emphasizing. At T = the decay width, T, 
depends on all the scalars in the problem. These are the masses M and m. 
The momenta of the decay photons, k\ and k 2 (ki + k 2 = q, the momentum 
of A) do not contribute anything else, since k\ ■ k 2 = M 2 . There are major 
changes in the kinematics at finite temperature, since the heat bath selects 
out a preferred frame. A covariant description can be retained in finite 
temperature field theory, provided an extra vector is introduced into the 
problem — the velocity of the heat-bath with respect to the frame in which 
one chooses to work, u (u 2 — 1) 0. The amplitude may now depend on all 
the scalars in the problem. There are four new scalars — the temperature 
T, the Fermion chemical potential /x, and two scalar-products, which can 
be chosen as q ■ u and k\ • u. Since T is obtained after integration over the 
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two- photon phase space, the dependence on k\ ■ u drops out. Consequently, 
T is a function of M, m, T, /i and q ■ u. This last quantity is the energy of 
A in the rest frame of the heat-bath. 

A simple representation of the new invariants facilitates the discussions 
in the later sections. Evaluate them in the rest frame of A, q — (M, 0, 0, 0). 
Align u along the ^-direction and write u = (7, 0,0, /3). Since the energy of 
each photon is M/2, and the momenta balance, 



(1-1) 



q ■ u = 


M 7 , 


k± ■ u — 


|M 7 (l-/3cosx), 


k 2 - u = 


iM 7 (l + /?cosx). 



It is now obvious that 7 and x are Lorentz invariants. Also, (3 = J (I — I/7 2 ) 
can be interpreted as the boost between the rest frames of A and the heat- 
bath. 

The plan of this paper is the following. In Section |2] we present the calcula- 
tion of the triangle diagram at T > 0, compute the decay width, and discuss 
some applications of these results. The next section contains a discussion 
of other related processes, including the Bose enhancement of multi-photon 
decays. Some field-theoretic points are discussed in the appendix. 

2 The Triangle Diagram at T > 0. 

In this section we compute the triangle diagram at T > 0. A limiting case, 
7 = 1, has been considered before 0. In all physical applications it is 
appropriate to consider the limit M <C m. The photon-Fermion coupling is 
taken to be iei(k) (for any vector p we denote p = 7^^), where e{k) denotes 
the polarisation vector of a photon of momentum k. The Fermion-A coupling 
is igj 5 . 

We summarise the computation of the decay width at T = 0. The matrix 
element, M., coming from the triangle diagram, has the form 

M = -^-teffay^T^fiM,™,) where T^ = me, uap k{k p 2 . (2.1) 

The Lorentz-scalar, /, is a form-factor for the effective A77 vertex and is 
obtained as an integral over the Fermion-loop momentum. It is a function 



only of M and m. It is possible to expand / as 

f(M,m) = f ( m )+z 2 f 1 (m) + ---, (*=|^)- (2-2) 

A text-book calculation gives /o = —7i 2 /m 2 and /i = — it 2 /3m 2 |J. The 
decay width, T, is computed by integrating the squared matrix element over 
the Lorentz-invariant phase space of the final state — 

r(M,m) = (^)- 2 M. (2.3) 

In the usual formulation of finite-temperature field theory , the physical 
fields (called type 1) are doubled by the addition of so-called "thermal ghosts" 
(type 2 fields). Each vertex involves only one type of fields. External legs 
connect only to type 1 vertices. Each propagator may connect any two types 
of vertices and hence becomes a 2 x 2 matrix which can easily be written 
down in terms of advanced and retarded Green's functions along with the 
Bose or Fermi distributions. In the triangle diagram all vertices connect to 
external legs, and we need to consider only the type 11 Fermion propagators 



T>(p,u) = (p + m) 



-n K + 2ni5(p 2 — m 2 ) x 

V -m +ie / 2 .4) 

{Q(po)F+(p ■ u) + Q(- Po )F-(p • «)} 



where a regulator has been placed in the denominator of the propagator, 
F ± (x) = l/(exp(|a;| =F y)/T + 1), and /x is the Fermion chemical potential. 
Note that this is the sum of a T = and a T > part. It is convenient to 
use a diagrammatic notation where the T = part is denoted by a line and 
the T > part by a line with a cut. Thermal effects on external legs can 
be subsumed into Bose and Fermi distributions multiplying the phase space 
volume element. We return to this point in the appendix. 

2.1 The Form Factor 

Note that the Dirac structure of the propagator is the same as at T = 0. As 
a result, the tensor structure of the matrix element remains as in eq. (|2[ 



and the only thermal effects come from the form factor / and the photon 
phase space. The modified form factor is 



/ = jd^pVi^Vip-k^Vip-q). 



(2.5) 



Separating the T = and T > parts of the propagator, the thermal form 
factor has eight terms. Apart from the T = contribution, other terms can 
be grouped as (a) one term with the thermal part of all three propagators 
(b) three terms with the thermal part of any two and (c) three terms with 
the thermal part of only one propagator. For M < 2m and m > 0, the three 
mass-shell conditions in (a) leave no phase-space volume to the integral and 
hence this contribution is identically zero. For the same reason each term 
in (b) is also zero. Consequently, the thermal contribution comes only from 
(c), written in a diagrammatic notation as 



./ 




+ q 




fJ>,ki\ 



v,k? 



(2.6) 



= 2(J! + J 2 + J 3 ). 

The factor of two comes from the diagrams obtained by simultaneous inter- 
change of (jl, v and k\, k%. 

Each of the integrals Jj can be reduced to a Lorentz invariant one- 
dimensional Fermi integral. We show some of the details for J\, to demon- 
strate that the computation can be performed covariantly After shifting 
p — > p + ki, and introducing a resolution of identity as a sum over positive 
and negative energy B-function, J\ reduces to a simple integral over a real 
particle phase space 



h 



-7T 



J 2p y 



x 



F + (x)) 



p- k\p ■ k 2 



(2.7) 



This integral is finite and well-defined for m > 0, and, being a Lorentz- 
scalar, may be evaluated in any convenient frame. We use the frame where 
u = (1,0,0,0). In this frame po = p ■ u is a Lorentz-scalar. As a result, 



p^ 



vl 



m 



is also invariant, and all the manipulations shown here are 



explicitly Lorentz invariant. 

First we use the Feynman trick to write 



h 



dpp 2 

Po 



F-(po)+F + (p )}n ll where Vt ± 



dadTt 



(2.8) 



Here a is the Feynman parameter (0 < a < 1) and V = aki + (1 — a)/c 2 . 
We choose the orientation of the frame by setting V = (Vo,0,0, V3). Since 
Vq = V ■ u, both Vo and V3 are Lorentz- invariant. The angular integral is 
easily performed, leaving Qi = (An/M 2 ) f da/Qi(a), where the Qi(a) is a 
quadratic in a with Lorentz-invariant coefficients. It is easy to check that 
the two roots of Qi do not lie in the range of integration. Then the integral 
over a is easily performed. With the notation 



V 2 



P 



2 1 2 2/i n2 

+ m 7 (1 — p 



cos' 



x), 



we have the compact expression 
47T 2 f°° dp p r 



Jx 



M 2 



PoV 



F-(p )+F + (p ) log 



'V + P 



(2.9) 



(2.10) 



J 2 and J3 can be obtained in a similiar fashion. 

The T > result may be exhibited in a form similiar to that in eq. (|2.2j) 
by expanding the integrals in a series in z = M/2m. Each of the integrals Jj 
(i = 1,2,3) has an overall factor of 1/M 2 . The only M dependence of J\ is in 
this factor. The expansion of J 2 + J3, starts at order 1/z. A straightforward 
check shows that the first term of the expansion precisely cancels J\. The 
first non-zero term is of order z° and this is followed by a series in z 2 . All 
these terms in the expansion come entirely from J 2 + J3. 

We give here the final result for the full T-dependent f — 



/ (M,m,/i,T,7) 
x 



7T 

m 



,22 



+ 2 HrJS°^ S [F-(Cro)+F+(Cro)] 



7r 7 roo drr 
(2r 2 - 7 2 /? 2 sin 2 X )\(l + P 2 ) log (g±T ) 



2rR 

r 2 + 7 2 



**$$ 



(2.11) 



(2r 2 /? 2 -(r 2 + 7 2 )/? 2 sin 2 x) 



where J3 = j3 cosx, r o = 1 + r2 > R 2 = r 2 + 7 2 (1 — J3 2 ), and C = m/T. This 
complicated looking result may be simplified in various limits, or evaluated 
numerically (see Figure [I]). 

An instructive limit is 7 — > 1, i. e., the special case when A is at rest in 
the rest-frame of the heat bath, and \x = — 



•2 



f(M,m,T) = - — + — — log J U— - -. 2.12 

m m Jo Tq \r — r/ e^ ° + 1 

This result agrees with previous expressions obtained in this limit ||. 

It is also interesting to see that in the limit T ^> m, the integral takes on 
the value 1/2 and the thermal part of / precisely cancels the T = result, 
giving / = 0. Thus the decay width vanishes in the joint limits (5 — > and 
M <^ m ^ T. The reason is quite simple. A close perusal of the T = 
calculation shows that / picks up a non-vanishing contribution only from 
that part of the phase space where one of the Fermions is on-shell. This 
contribution is fully Pauli-blocked in the limit T ^> m. 

For m <C T it is useful to perform an expansion of the form factor /. 
From eq. (|2.12|) it is clear that 



0rf iog (^) +oK2) -°- 77o7c S 



In the other limit, m ^> T, the thermal contribution vanishes exponentially 
in C, leaving / to take on its T = value. 

The results are even simpler when A travels very fast through the medium, 
i.e., in the limit 7 — ► 00. In this limit the quadratics appearing in the 
definitions of all the J's go as 7 2 . Hence the thermal contribution vanishes 
as I/7 2 , and the matrix element is equal to its T = limit. 

2.2 The decay width 

Next we construct the decay width for A — > 77. In addition to the matrix 
element computed in the previous section, we require the 2-body decay phase 
space for T > 0. This contains a factor 1 + B(k ■ u) for each photon (due to 
stimulated emission in the heat-bath), where B(x) = l/(exp |x| — 1) is the 
Bose distribution function. Then, after the usual reduction of the 2-body 
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Figure 1: The form-factor \f Q m 2 /7i 2 \ shown as a function of 7 and cosx for 
values of T and il appropriate to a supernova core (T = 40 MeV, 11 = 450 
MeV; lower panel) and the neutrino-sphere in a supernova (T = 1 MeV, 
fx = 0; upper panel). Isolines drawn on the 7-cosx plane correspond to 
function values increasing in steps of 0.05 upto a largest value (on the right) 
of 0.95. 



phase space, 



III J\/fO 

r(M, m, //, T, 7) = 9 j^~7 — I-i d cos xf 2 (M,m,fj,,T, 7, x) 

[l+B(A;i-«)][l+S(A;2-«)]. 

(2.14) 
For the value of / calculated in the previous section, the integral has to be 
performed numerically for general values of the argument. 

When the x~ dependence of the form factor can be neglected, the thermal 
effect arises entirely from the phase space factors. In this case the decay 
width is 



„,„, m , g 2 a 2 m 2 M 3 f 

Y{M, m, //, T, 7) = 



64/?^7T 7 Vl-e- 2 « 



/ e (i+/3)S _-T 
/3£-log 



eC 1 "^ - 1 

(2.15) 
where £ = M7/2T. For /? -> and M < T, we find that T is the T = 
width times T 2 /M 2 . On the other hand, for TcMwe recover the T = 
result. For general values of the parameters a numerical integration of the 
expression in eq. ( |2.14j ) is required. 



The physics of the temperature dependence of the decay width can be 
summarised as — 

1. For T C If < m or large (3, thermal effects can be neglected and the 
decay width is the same as at T = 0. 

2. For small (3 and M <^ T ^C m the dominant effect is the stimulated 
emission of final state photons, leading to a highly enhanced decay rate. 

3. For small (3 M <C m <C T Pauli blocking of the on-shell part of the 
loop determines the physics and the decay rate is reduced. 

2.3 Some Applications 

The simple computations above can be used to bound M. A large effect 
comes from the fact that V ~ M 3 at T = 0, whereas V ~ T 2 M for M < T. 
This thermal enhancement of the decay width has an interesting consequence 
when M -C Tfe = 0.235 x 10~ 3 eV. The T = result, in conjunction with 
some cosmological arguments, would predict a bound on g < c/M 3 ^ 2 , whereas 



the T > results strengthen this to g < d /v M (c and d are computable 
constants). 

If A is to be a candidate for cold dark matter, then (3 « and the lifetime 
must be greater than the age of the universe, to [@- Then 



where the average Fermion mass m~ 2 = J2f m l 2 is dominated by the 
lightest flavour, the electron. Assuming the present age of the universe to be 
15 x 10 9 years, we find that when M 3> T b , c = 1.2 x 10~ 7 eV 3 ^ 2 and when 
M < T b , d = 2.5 x 10" 4 eV 1/2 . 

The coupling g also can be constrained by the condition that the energy 
density of the decay photons is limited by observation of the micro-wave 
background, p 7 . We will assume || that A is produced non-thermally in the 
early universe. Then we can use the decay width in the limit /3 — *■ 0. The 
additional photon energy density due to the decay of A at present is 

Ap 7 (t) = ~p A {t) [l - e- r ^>] , (2.17) 

where ti is the production epoch and p A is the energy density which would 
have remained with A in the absence of decays. Now Ap 7 < p 7 , and ~p A must 
be greater than the present critical density p cr M. Since ti <^C t, we find 



g < (—) (l - e-W) HWl-fh/tPcr) (218) 

Taking p A = zp cr with z = 0.9, p 7 = 0.260 eV/cm 3 , p cr = 1.054 x 10~ 5 /iq 
GeV/cm 3 and choosing /i = 0.85 (to get the weakest bound), we find that 
c = 7.1 x 10~ 10 eV 3/2 and d = 1.5 x 10" 6 eV 1/2 . 

On the other hand, if we assume that Ap 7 at the present epoch is smaller 
than the known error induced through the error on the COBE measurement 
of Tb, then the constraints are harder by a factor of 0.04. This gives c = 
2.8 x 10~ n eV 3/2 and d = 6 x 10~ 8 eV 1/2 . In all these cases the bound 
on g is sharpened by two orders of magnitude through the leading order 
diagram alone. Although this does not place very stringent restrictions on 
axion models, in the next section we consider higher order diagrams and 
argue that they may lead to even sharper bounds. 
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Apart from axion models, these computations can also be used to bound 
the couplings and masses of any other pseudo-scalar with a 75 vertex to 
charged Fermions. One such model, which is not ruled out by data, is a 



singlet-triplet Majoron model [K|. This Majoron is a Goldstone Boson for 
Lepton number breaking and may be given a mass through couplings to 
gravity ]TT[ . Since this mechanism typically produces low masses, the thermal 
bound is applicable. 

Another class of limits arises from consideration of the cooling rates of 
stars from emission of pseudo-scalars. Our computation of the form factor is 
applicable to this situation, not only through the direct two-photon channel, 
but also through the Primakoff process. Thermal effects in most stars are 
rather small, due to the fact that T <C m. However, they can be important 
in supernovae since the A77 vertex is suppressed (see Figure |l|). This may 
affect bounds on axions from supernova cooling rates|rj, |T3] ]. 

3 Other Processes. 

3.1 Multi-photon decays 

We have computed thermal effects on the decay width of a pseudo-scalar 
(of mass M) in the two-photon channel through the triangle diagram with 
charged Fermions of mass m circulating in the loop. We found a strong 
increase in the decay width due to Bose enhancement of the final state pho- 
tons for M < T <ti m, leading to a sharpening of the bound on the cou- 
pling between the pseudo-scalar and Fermions by two orders of magnitude 
for M = 10~ 5 eV. Multi-photon decays are likely to yield more stringent 
bounds. 

The total width of the axion into multi-photon channels at T = can be 
written as the perturbation series 

T = a 2 T 2 + a 4 r 4 + • ■ ■ (3.1) 

Since a 2 ~ 10~ 4 , the series converges extremely fast and even the second 
term can be neglected. However, the situation changes at finite temperature. 
The higher order terms, involving multi-photon decays are enhanced due to 
Bose factors on each external photon leg, and some of these may compensate 
the damping due to the fine structure constant. 
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Examine the case of a pseudo-scalar decaying into In photons, where the 

i-th photon momentum is hi. When each ki ~ M/2n, the Bose factors go as 

(2nT/M) 2n . As a result, the contribution of this process to the total decay 

width is 

r /2nT\ 2n /M\ 2n 



a 2n T 2n = g 2 a 2n -^- — M — ■ ( 3 - 2 ) 

y (2n)\ \ M ) Km) K J 

The last factor comes from propagators and the Dirac trace, the factorial is 
due to photon counting and c n is essentially the dimensionless contribution 
due to the Fermion loop integral. Clearly, this region of phase space is not 
very important at higher orders since T Cm. 

However, in the part of phase space where one or more ki — > the Bose 
factors can become arbitrarily large. This can happen for n > 2. It is easy 
to check that the real 2n photon emission diagrams give non-vanishing am- 
plitudes in this region of phase space. This is not the usual Bloch-Nordsieck 
problem, since the T = width is perfectly well defined and insensitive to the 
infra-red (as long as m > 0). We expect these putative divergences to cancel 
when thermal ghosts and virtual corrections are taken into account. The 
result is that, the decay width at T > is infra-red sensitive and hence must 
be resummed over all numbers of photons. Demonstrating the cancellation 
and extracting the finite parts, prior to summing the series then requires the 
full power of thermal field theory. This work is under progress and will be 
reported elsewhere. 

3.2 Stellar Cooling 

Stellar cooling arguments are dependent on the coupling of axions to matter. 
For small couplings axions stream out of the star, carrying energy. At larger 
couplings, the mean-free path may be smaller than the radius of the star, 
and an axio-sphere may form. Axions escape only from the surface of the 
axio-sphere. 

We have seen that the axion-photon form factor decreases significantly at 
temperatures appropriate to the core of supernovae. This does not have a sig- 
nificant effect on the cooling rate due to low-mass axions, since the decrease 
affects only the low-energy end of the spectrum. For small M the energy loss 
due to low-energy axions is negligible and large changes in production rate 
can be tolerated. 
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The main effect is seen for M large enough that an axio-sphere may be 
formed. Apart from the decrease in the form factor, one should also take into 
account changes in the pion decay constant, /„., and the pion mass, m f , since 
they enter the axion coupling. At T « 40 MeV, chiral perturbation theory 
shows some decrease in both these quantities [|14J]. Further, suppression of 



interaction rates of low-energy axions with leptons or hadrons, due to Pauli 
blocking, increases the mean-free path of low-energy axions. The net effect 
is that low-energy axions may escape from the core. At the same time, high- 
energy axions are converted to low energies by processes such as Ae — > Ae 
(Compton). The inverse process of boosting the energy of the axions is 
suppressed by Pauli blocking and the falling density of low-energy axions. 
This may lead to a destabilisation of the axio-sphere. Quantitative estimates 
of the relaxation time for this instability and the equilibrium phase space 
distribution of axions will be detailed elsewhere. 

Even in parts of parameter space where this instability is negligible, es- 
timates of the axio-sphere radius and the critical coupling at which it forms 
necessitate a full thermal field theory computation. The radius of the axio- 
sphere is closely related to the relaxation time (and the mean free path 
length) of the axion. The relaxation rate is the difference of the production 
and decay rates (see the appendix) and is given by the imaginary part of the 
axion two-point function. This quantity is strongly affected by hard thermal 
loops, if the mass of the axion is less than yfaT '. Assuming T « 40 MeV, 
treatments of axions with mass less than about 4 MeV in the supernova core 
require hard thermal loop resummation. 

A Some Field-theoretic Niceties. 

In Section |2| we have written down rules for computation of T > decay 
widths which look very similiar to those at T = 0. The formalism for doing 
this is developed in |L5], [Tj|. We summarise the results in this appendix and 
show that our techniques are justified. 

In general one computes decay rates by writing down cutting rules for 
loop contributions to a two-point function. At T > loop diagrams require 
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the matrix propagator JFJ 

iG{p) = U(T,p)(f^ s ? (p) )u(T,p), 
S(p) - 



i 

72 ZT 



p -m +ie (A.l) 



tG ± (p) = 2n5(p 2 -m 2 )[Q(±p ) + B(p-u)]. 

These rules for Bosons can be generalised in a similiar form to Fermions and 
Gauge Bosons. At T = the matrix is diagonal and only the G(±po) term 
appears on G ± . 

The T = notion of cutting a line is generalised at T > to the notion 



of "circling" a vertex |15|]. If a G n line is circled only at one end, then it is 
replaced by G + if the momentum flows to the circled vertex (G~ otherwise). 
For a G22 line the rule is reversed, and the off-diagonal lines are left un- 
touched. Since the off-diagonal terms are absent at T = 0, only 0-functions 
appear on such lines, and these reproduce the familiar cuts. 



Furthermore [lq|, the imaginary part of the self-energy can be written as 



Imll(p) = p-uT r (p-u) = p-u(T — T p ), (A. 2) 

where T r is the relaxation rate, T p is the production rate of the particle in 
the heat bath and T is the computed decay rate. 

Now, for the decay rate, we require that the particles in the loop are 
observed in the final state. We might guess that the computation requires 
the replacement 

iG ± (p) = 2TT5(p 2 -m 2 )&(±p )[l + B(p-u)}. (A.3) 

The overall factor of G allows us to talk of cut lines as at T = 0. Similarly, 
for the production rate we might guess that the replacement, 

iG ± (p) = 27r5(p 2 -m 2 )Q(TPo)[l + B(p-u)}, (A.4) 

is called for. Then do we miss cross terms which involve a 8(po)8(—p' )? A 
simple check reveals that if the two-point function involves a real on-shell 
particle, then such cross terms vanish due to kinematical reasons. This is the 
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content of eq. ( |A.2| ), and the justification for the prescription in eq. ( |A.3| ), 
used in Section |[ 

In our computation of the decay width, T, we have neglected hard thermal 
loops ||, corresponding to screening of electric charges. Knowing that they 
are important for computations of the photon self-energy at T > 0, for 
k ■ u ~ eT, is their neglect justified? The answer is: yes, because these have 
large contributions to T r , but their effect on T p and T separately are small. 
We give the outline of a demonstration below. 

The hard thermal loop contributions to the process we are interested in 
can be represented as 

p(a6) _ J a ^ L (A. 5) 

The blobs represent all vertex corrections to the triangle diagram and self 
energy corrections to the Fermion and photon lines. The result is, retaining 
only the terms of interest, 

r (a6) = J d 4 hd 4 k 2 M^M^ lmU^(k x ) lmU vp (k 2 ) 

~ / d A k x d A k 2 [1 + B{k x ■ u)\ [1 + B(k 2 ■ u)} 

x M^M$ ) p l {k 1 )p l {k 2 ) + M% ) Mi?p l {k 1 ) Pt {k 2 ) 

+M^M$ ) p t {k 1 ) Pl {k 2 ) + M^M^p^kJptih 



(A.6) 



where pi and p t are the spectral densities of the longitudinal and transverse 
components of the photon. The quantities M^'^ are the contributions from 
the blobs; the subscripts denote appropriate components of the tensors. The 
analysis of |I7| can be adapted to show that the term quadratic in p t is the 



dominant infra-red term and gives the same contribution as that obtained 
with tree level propagators for the photon. Hence, hard thermal loop re- 
summation does not change our results for T. Similiar arguments hold for 
the production rate V p . However, the leading terms, quadratic in p t , cancel 
between the production and decay rates and hence the sub-leading terms 
become important for a reliable computation of the relaxation rate r. r . For 
this reason, hard thermal loop resummation becomes important for the com- 
putation of the radius and other properties of the axio-sphere. 
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